Dynamics of optically generated vortices in a one— component ultracold fermionic gases 
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We show that the phase imprinting method is capable of generating vortices in a one-component 
gas of neutral fermionic atoms at zero and finite temperatures. We find qualitative differences 
in dynamics of vortices in comparison with the case of Bose-Einstein condensate. The results of 
the imprinting strongly depend on the geometry of the trap, e.g., in asymmetric traps no single 
vortex state exists. These observations could be considered as a signature of Cooper-pair based 
superfluidity in a Fermi gas. 
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One of the spectacular properties of a macroscopic 
bosonic-type quantum system is a superfluidity, origi- 
nally observed in a liquid helium II. Although superflu- 
idity is a complex phenomenon, it is certainly intimately 
related to the existence of quantized vortices. The recent 
experimental realization of Bose-Einstein condensation 
in confined alkalai-metal gases allowed to study this 
connection in details. Quantized vortices (in a form of 
small arrays as well as completely ordered Abrikosov lat- 
tices) have been already observed in many laboratories 
[PJ. Interferometric detection method showed directly 
2-7r-phase winding associated with the presence of a vor- 
tex in the condensate |||. Of course, this is a manifesta- 
tion of the macroscopic wave function. 

However, the degenerate Fermi gas has no macro- 
scopic wave function. It is rather described (at zero 
temperature) in terms of the many-body wave function, 
built of single-particle orbitals, whose evolution is gov- 
erned by the Schrodinger equation. Hence, the appear- 
ance of quantized vortices in a fermionic system demon- 
strates the BCS-type transition from normal to super- 
fluid phase. Recently several groups have undertaken the 
effort to achieve quantum degeneracy in a dilute Fermi 
gas In JILA experiment j| §] 40 K atoms were 

trapped in two hyperfine states and cooled evaporatively 
by collisions between atoms in a different spin states 
whereas in Refs. || ^| a mixture of bosonic ( 7 Li) and 
fermionic ( 6 Li) atoms was used to perform the sympa- 
thetic cooling of fermions. 



In this Letter we investigate the dynamics of optically 
generated vortices in a Fermi gas in a normal phase and 
compare it with the properties of vortices in the Bose- 
Einstein condensate. A good approximation to the spin- 
polarized Fermi gas at low temperatures is that of non- 
interacting particles (the s-wave scattering is absent for 
spin-polarized fermions). We then assume that at zero 
temperature the system is described by the Slater de- 
terminant with the lowest orbitals being occupied. The 
details of the technique of phase imprinting are thor- 
oughly discussed in Refs. [|| || for bosons and repeated in 
Refs. [|l0, 11 for fermions. It is worth stressing that this 
method was already experimentally realized for the gen- 
eration of solitons in Bose-Einstein condensate jllL |l3| . 

We start our considerations with a two-dimensional 
symmetric trap and imprint the phase at the center of the 
trap. Each atom acquires a phase given by k(j), where k 
is an integer number and cj> is the azimuthal angle around 
the origin (each wave function is multiplied by the factor 
exp(jfc(/>)). The following free evolution of the single- 
particle orbital can be calculated by using the propagator 
technique 



ipj(x,y,t) 



K(x,y,x',y',t)ifj(x',y',0)dx'dy' 



(1) 

where K(x,y,x',y',t) is the propagator function for 
two-dimensional symmetric harmonic oscillator pM 



K(x,y,x',y',t) = ——. exp <^ — [{x 2 + y 2 + x' 2 + y' 2 coswi- 2xx' - 2yy'] \ . 
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In this case it is convenient to represent the 
single-particle orbitals in polar coordinates (hereafter, 
y/h/ (Mu) and l/u> are used as units of length and time, 

respectively): Vnm (r, 0) = A nm rl m l l)™\r 2 ) e~ r ' ''/ 2 e™"> 
and A nm = (n\/[n(n + |m|)!]) 1//2 and the integration ac- 



cording to the prescription (|lj) can be performed first over 
the azimuthal angle then over the radial distance. 
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The first stage requires the use of the formula pi 



-iz cos @ 



dO = (-i) l 2<K Ji{z) , 



(2) 



where Ji(z) is the Bessel function of integer order and z 
any complex number. In our case z is real and equals 
rr'/ s'mt. In the second stage one has to calculate 



J m+k (br')L\f{r' 2 )dr\ (3) 



where a 2 = 1/2 — i cost/ (2 sini) and b = rj smt. Since 
the generalized Laguerre function L'n is a polynomial, 
the expression (^) can be expanded and integrated term 
by term with the help of appropriate formula for integrals 
of Bcsscl functions Ml 



./<M-l) e -aV j (br > )dr > = ^±MM)1 



xAf [-,+ -/,, + 1,- — 



(4) 



valid for Re(a 2 ) > and Re(/i + z/) > 0. M(a,j,z) is a 
confluent hypergeometric function. 

The result is the product of b m+k and the combination 
of the confluent hypergeometric functions; each of the 
argument equal to — 6 2 /(4a 2 ) 

A (_A\rn+k-l 

(p(r, <(>,t) = — e 



sint 



T(l + m + k) 2 1 + m + k 
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(-1) ( n + \m\ U(l + ! + 



1=0 

x M | 1 + I 



n — I 



m+/c+|m| 



1 2+2l+m+k+\m\ 



m + k + \m\ , b 2 
2 ±A > l + m + k >-^ 2 



Hence, the value of the wave function (single-particle or- 
bital) at the origin is fully determined by the factor b m+k 
and equals zero provided the phase imprinting is strong 
enough (k is larger than the maximum absolute value of 
magnetic number m of occupied orbitals). In such a case 
the diagonal part of the one-particle density matrix is 
zero at the origin all the time. On the other hand, the 
part of the phase of particular orbital dependent on the 
azimuthal angle is given by (m + k)<p and the circulation 
along the circle centered at the vortex line (point) is 



M 



k) 



£ 7 >i(r,t)| s 



In Fig. [j] we plotted the circulation for various paths 
(circles centered at the vortex point) at different times 



for a system of N — 6 fermions after writing the phase 
3(f>. It is clear that the strong enough phase imprinting 
leads to the excitation of the Fermi gas with non-zero 
(but not quantized) circulation around the point where 
the density vanishes. We call this state the vortex. 
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FIG. 1: Circulation along the circle as a function of its radius 
for various moments: 0.05 (solid line), 0.26 (dashed line), and 
1.57 (dotted line) in units of 1/uj. The system and phase 
imprinting parameters are as follows: TV = 6 and k — 3. 



For weaker phase imprinting (e.g., for k = 2 and 
N = 6) the response of the system becomes more com- 
plex. The density just does not go to zero at the center. 
In some cases the density at the origin oscillates between 
zero and local maximum and the circulation goes to zero 
while the radius of the circle decreases (vortex-antivortex 
oscillations). In others, the dip in the density is strongly 
diminished, however the vortex reappears at multiple in- 
tegers of the trap period. 

It turns out that we can follow analytically the evolu- 
tion of the single-particle orbitals also in the case when 
the absorption plate (and the laser beam) is shifted off 
the center of the trap. For that, the single-particle or- 
bitals arc multiplied by the factor 



(x - s ) + i(y - 2/q ) 
y/(x- x ) 2 + (y- y ) 2 



where (xo,yo) are coordinates of the laser beam with re- 
spect to the center of the trap. This time it is better to 
represent the atomic wave functions in Cartesian coordi- 
nates: ip nin2 (x,y) = A ni A n2 H ni (x)H n2 (y)e-( x2+y2 ^ 2 , 
where A n = (l/( 7 r 1 / 2 2 Tl n!)) 1 / 2 and H n (x) is the n-th 
Hermite polynomial. 

To find the evolution of single-particle orbitals, gov- 
erned by Eq. (Q), it is convenient to shift the origin of 
Cartesian coordinates (x',y') by the vector (xo,j/o) an d 
then introduce the polar coordinates. The integration 
over the azimuthal angle requires the use of generating 
function technique and the following extension of formula 
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(v^ 2 + io 2 ) fc 



with z and w being complex numbers. The result 
(Inin 2 ( r ')) f° r an y state (nx,*n.2) is obtained by expand- 
ing the left-hand side of the equation 



oo oo 



(5) 



where 



m— n,2— 



z = — — (zoe 1 * — x) — i2si 
sini 

w = (y Q e lt - y) - i2s 2 

smt 



in a Taylor series with respect to (si,S2) around (0,0) 
point. 
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FIG. 2: Evolution of the density distribution of N = 105 
atoms in two-dimensional symmetric harmonic trap after im- 
printing a phase of (14 0) 1.5 osc. units off the center. The 
successive frames correspond approximately to the moments: 
(a) 0, (b) 1/4, (c) 1/2, and (d) 3/4 of trap period. The vortex 
rotates around the center of the trap changing periodically 
the size of its core. 

Since the only dependence on the distance r' appears 
through the J Bessel function it is convenient to perform 
the integration over the radial distance first. Again, by 
using the formula (|J) one gets (up to time-dependent 
constant) 



for the left-hand side of (||). The expression (^) implies 
that, assuming the phase imprinting is strong enough 
(k > max(ni,7i2)), each atomic wave function includes 
the factor (z + iw) in an integer power higher than or 
equals 1, taken at the point (si = 0, s 2 = 0). It is easy 
to check that 

(z + iw) Sl=0t s 2= o = (x cost - j/o sini - x) + 

i (xq sin t + yo cos t — y) . (7) 

It means that for k > max(ni,n2) all orbitals have zero 
at the point (z + iw) Sl= o,s 2 =o = 0. So, this zero is also 
present in the diagonal part of the one-particle density 
matrix. Moreover, it is clear from (0) that the zero ro- 
tates around the center of the trap at the radius equal 
to \/ x^ + y 2 with constant frequency which is the trap 
frequency. For n\ = n 2 = the time-dependent wave 
function is given by 



r(i + f) AqAqi*- 1 



r(l + fc) 2 1 + k a 2 + k sint 



=«(- 2 +^o)-^o 2 



4a 2 J ' 



where z\ = (xq exp (it) — x)j sini, W\ = (yoexp(it) — 
y)/ smt, and b 2 = z\ + w\. Evolution (obtained numer- 
icallyby solving the many-body Schrodinger equation 
Jl0| , pd| ) of the density distribution for N — 105 atoms is 
shown in Fig. |^. 
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-(sl+sl) e 2{ Sl x + S2 y ) ( z + iw )k 



(6) 
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FIG. 3: Density distribution (cuts along 'x' axis) of N — 105 
(average number) atoms at temperatures zero and O.ITf af- 
ter writing the phase 140. The curves correspond to moments 
0.75 (thick and thin solid lines for O.ITf and zero tempera- 
ture, respectively) and 1.5 (thick and thin dashed lines for 
O.ITf and zero temperature, respectively) in units of 1/w 
and are the results of averaging procedure over 1000 config- 
urations. The thick dashed curve shows the moment of the 
lowest contrast (which is about 50%). 

We consider now the generation and dynamics of vor- 
tices at finite temperatures. To this end, we employ the 
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Monte Carlo algorithm and generate a number of many- 
body configurations JTl| according to the Fermi-Dirac 
statistics. At finite temperatures the single-particle 
states with the energy above the Fermi level are pop- 
ulated. It may happen that for a given phase imprint k, 
even strong enough to generate a vortex at zero temper- 
ature, there are configurations with such single-orbitals 
that m + k = 0. For these configurations we will see 
again complex structures like vortex-antivortex oscilla- 
tions in the center of the trap or circulating around it, 
depending on the way the phase imprinting was done. 
Averaging procedure, according to the grand canonical 
ensemble rules just shows the vortex with the lower con- 
trast (see Fig. |J). However, one might argue that in 
a particular experimental realization only one configura- 
tion is involved and hence vortices, reappearing vortices 
or vortex-antivortex oscillations should be observed. In 
any case, non of these structures is dissipatively circulat- 
ing out of the system as predicted for the vortices in the 
Bose-Einstein condensate at finite temperatures jig. 
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FIG. 4: Evolution of the density distribution of N = 105 
atoms in two-dimensional asymmetric harmonic trap (30Hz x 
50Hz) after imprinting a phase of (14 0) at the center. The 
successive frames correspond approximately to the moments: 
(a) 0, (b) 1/4, (c) 1/2, and (d) 3/4 of trap period. The vortex 
reappears at integer multiples of trap period. 

Finally, we have done some numerical calculations for 
asymmetric traps at zero temperature. The results are 
presented in Fig. |]. Again, striking difference in com- 
parison with the Bose-Einstein condensate case is found. 
No vortex is generated, presumably because there is no 
common zero for orbitals after phase imprinting in the 
case of asymmetric trap. 

In conclusion, we have shown that the dynamics of 
vortices generated in a cold gas of neutral fermionic 
atoms significantly differs from the case when they are 
generated in the Bose-Einstein condensate. At zero 
temperature and symmetric harmonic trap, the vortex 
rotates around the center of the trap. However, for 
fermions it periodically changes the size of its core. For 
bosons, the size of the vortex core remains at the order 
of the healing length. Further differences appear at 
finite temperatures. Contrary to dissipative escape of 
vortex from the trap predicted for the Bose-Einstein 
condensate fl6}| , we find the stable rotational motion of 
the fermionic vortex. Finally, the case of asymmetric 
trap shows new density patterns present in a Fermi gas 
and not observed for bosons. All these findings could be 
used as a test of BCS type phase transition in a Fermi 
gas. 
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